CLNS 04/1901 
CLEO 04-19 



(N 



Improved Measurement of the Form Factors in the Decay 

A+ -»■ Ae+v P 



J. W. Hinson, 1 G. S. Huang, 1 J. Lee, 1 D. H. Miller, 1 V. Pavlunin, 1 R. Rangarajan, 1 
B. Sanghi, 1 E. I. Shibata, 1 I. P. J. Shipsey, 1 D. Cronin-Hennessy, 2 C. S. Park, 2 
W. Park, 2 J. B. Thayer, 2 E. H. Thorndike, 2 T. E. Coan, 3 Y. S. Gao, 3 F. Liu, 3 
R. Stroynowski, 3 M. Artuso, 4 C. Boulahouache, 4 S. Blusk, 4 E. Dambasuren, 4 
O. Dorjkhaidav, 4 R. Mountain, 4 H. Muramatsu, 4 R. Nandakumar, 4 T. Skwarnicki, 4 

S. Stone, 4 J.C. Wang, 4 S. E. Csorna, 5 I. Danko, 5 G. Bonvicini, 6 D. Cinabro, 6 
M. Dubrovin, 6 S. McGee, 6 A. Bornheim, 7 E. Lipeles, 7 S. P. Pappas, 7 A. Shapiro, 7 
W. M. Sun, 7 A. J. Weinstein, 7 R. A. Briere, 8 G. P. Chen, 8 T. Ferguson, 8 G. Tatishvili, 8 
H. Vogel, 8 M. E. Watkins, 8 N. E. Adam, 9 J. P. Alexander, 9 K. Berkelman, 9 V. Boisvert, 9 
D. G. Cassel, 9 J. E. Duboscq, 9 K. M. Ecklund, 9 R. Ehrlich, 9 R. S. Galik, 9 L. Gibbons, 9 



B. Gittelman, 9 S. W. Gray, 9 D. L. Hartill, 9 B. K. Heltsley, 9 L. Hsu, 9 C. D. Jones, 9 
J. Kandaswamy, 9 D. L. Kreinick, 9 A. Magerkurth, 9 H. Mahlke-Kriiger, 9 T. O. Meyer, 9 
N. B. Mistry, 9 J. R. Patterson, 9 D. Peterson, 9 J. Pivarski, 9 S. J. Richichi, 9 D. Riley, 9 
A. J. Sadoff, 9 H. Schwarthoff, 9 M. R. Shepherd, 9 J. G. Thayer, 9 D. Urner, 9 T. Wilksen, 9 
A. Warburton, 9 M. Weinberger, 9 S. B. Athar, 10 P. Avery, 10 L. Breva-Newell, 10 V. Potlia, 10 
H. Stoeck, 10 J. Yelton, 10 K. Benslama, 11 C. Cawlfield, 11 B. I. Eisenstein, 11 G. D. Gollin, 11 

I. Karliner, 11 N. Lowrey, 11 C. Plager, 11 C. Sedlack, 11 M. Selen, 11 J. J. Thaler, 11 
J. Williams, 11 K. W. Edwards, 12 D. Besson, 13 S. Anderson, 14 V. V. Frolov, 14 D. T. Gong, 14 
Y. Kubota, 14 S. Z. Li, 14 R. Poling, 14 A. Smith, 14 C. J. Stepaniak, 14 J. Urheim, 14 
Z. Metreveli, 15 K.K. Seth, 15 A. Tomaradze, 15 P. Zweber, 15 S. Ahmed, 16 M. S. Alam, 16 
J. Ernst, 16 L. Jian, 16 M. Saleem, 16 F. Wappler, 16 K. Arms, 17 E. Eckhart, 17 K. K. Gan, 17 
Q_i| C. Gwon, 17 K. Honscheid, 17 H. Kagan, 17 R. Kass, 17 T. K. Pedlar, 17 E. von Toerne, 17 

H. Severini, 18 P. Skubic, 18 S.A. Dytman, 19 J.A. Mueller, 19 S. Nam, 19 and V. Savinov 19 

(CLEO Collaboration) 

^ ! 1 Purdue University, West Lafayette, Indiana ^7907 

2 University of Rochester, Rochester, New York 14627 
3 Southern Methodist University, Dallas, Texas 75275 
^Syracuse University, Syracuse, New York 13244 
5 Vanderbilt University, Nashville, Tennessee 37235 
6 Wayne State University, Detroit, Michigan 48202 
7 California Institute of Technology, Pasadena, California 91125 
8 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 
9 Cornell University, Ithaca, New York 14853 
10 University of Florida, Gainesville, Florida 32611 
11 University of Illinois, Urbana- Champaign, Illinois 61801 
12 Carleton University, Ottawa, Ontario, Canada K1S 5B6 
and the Institute of Particle Physics, Canada 
13 University of Kansas, Lawrence, Kansas 66045 
^University of Minnesota, Minneapolis, Minnesota 55455 
15 Northwestern University, Evanston, Illinois 60208 
16 State University of New York at Albany, Albany, New York 12222 



1 



11 Ohio State University, Columbus, Ohio 43210 
18 University of Oklahoma, Norman, Oklahoma 73019 
19 University of Pittsburgh, Pittsburgh, Pennsylvania 15260 
(Dated: December 29, 2004) 

Abstract 

Using the CLEO detector at the Cornell Electron Storage Ring, we have studied the distribution 
of kinematic variables in the decay — > Ae + f e . By performing a four-dimensional maximum 
likelihood fit, we determine the form factor ratio, R = /2//1 = —0.31 ± 0.05(stat) ± 0.04(syst), the 
pole mass, M po i e = (2.21 ± 0.08(stat) ± 0.14(syst)) GeV/c 2 , and the decay asymmetry parameter 
of the A+, a Ac = -0.86 ± 0.03(stat) ± 0.02(syst), for (q 2 ) = 0.67 (GeV/c 2 ) 2 . We compare the 
angular distributions of the A+ and A c and find no evidence for CP-violation: Aa c = |°^ c ^^ c j = 

0.00 ± 0.03(stat) ± 0.01 (syst) ± 0.02, where the third error is from the uncertainty in the world 
average of the CP- violating parameter, Aa, for A — > pir~ . 
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The charm quark is unstable and decays via a first order weak interaction. In semileptonic 
decays, which are analogs of neutron (3 decay, the charm quark disintegrates predominantly 
into a strange quark, a positron and a neutrino. The rate depends on the weak quark 
mixing Cabibbo-Kobayashi-Maskawa (CKM) matrix element |V^ S | and strong interaction 
effects, parameterized by form factors, which come into play because the charm quark is 
bound with light quarks to form a meson or baryon. Charm semileptonic decays allow a 
measurement of the form factors because |V^ S | is tightly constrained by the unitarity of the 
CKM matrix Q. 

Within the framework of Heavy Quark Effective Theory (HQET) j2j, semileptonic (J p = 
l/2 + — ► l/2 + ) transitions of A-type baryons are simpler than mesons as they consist of a 
heavy quark and a spin and isospin zero light diquark. This simplicity leads to more reliable 
predictions |3|, |4j for form factors in heavy-to-light transitions. The measurement of form 
factors in the A+ — > Ae + u e transition provides a test of HQET predictions in the charm 
baryon sector, a test of Lattice QCD, and information for the determination of the CKM 
matrix elements \V c b\ and \V u b\ using A° decays since HQET relates the form factors in A+ 
semileptonic decays to those governing A° semileptonic transitions. 

In the limit of negligible lepton mass, the semileptonic (J p = l/2 + — > l/2 + ) transition 
of a A-type baryon is parameterized in terms of four form factors: two axial form factors, 
Ff and F 2 A , and two vector form factors, FY and F%- These form factors are functions of 
q 2 , the invariant mass squared of the virtual W + . The decay can be described in terms of 
helicity amplitudes H\ A \ W = H\ A \ W + H XaX , where Aa and \\y are the helicities of the A 
and W + . The helicity amplitudes are related to the form factors in the following way 

lq 2 HY, = yjQ- [(M Ac + M A )FY - q 2 F 2 h 




Hi = J2Q_ [-FY + (M Ac + M A )if ], (1) 



!q*Hf Q = [(M Ac - M A )F* + q 2 F 2 % 

A I -Mr \ , , 77V1 



m x = pQ + [-Ft - (M Ac - M A )F 2 V ], 
where Q± = (M Ac ± M A ) 2 — q 2 . The remaining helicity amplitudes can be obtained using 

angular distribution can be written as [4] 5] 



the parity relations HY^Lx — +( — )^al^w In terms of the helicity amplitudes, the decay 



"n YD { A — N ^ | T r |2 q 2 P 

S = dq 2 d cos 9 A d cos 9 w d X = ( ~* W '2 (2^' cs| 24M^ * 

{f (1 - cos 9 W ) 2 \H il | 2 (l + a A cos# A ) + -(1 + cos9 w ) 2 \H_i_ 1 \ 2 (l - a A cos# A ) (2) 

o 2 O 2 

3 

+ - sin 2 ^vf[|-^|o| 2 (1 + a A cos# A ) + |if_i | 2 (l — a A cos# A )] 
3 

y=a A cos^sin^vK sin# A [(l — cos 9w)Re(H~ lo-^ii) + (1 + cos 9w)Re(Hi H^ 1 

2\/2 22 2 2 

where Gf is the Fermi coupling constant, \V CS \ is a CKM matrix element, P is the magnitude 
of the A momentum in the A+ rest frame, 9 A is the angle between the momentum vector of 
the proton in the A rest frame and the A momentum in the A+ rest frame, 9w is the angle 
between the momentum vector of the positron in the W + rest frame and the A momentum 
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in the A+ rest frame, x is the angle between the decay planes of the A and W + , and a a is 
the A — > decay asymmetry parameter measured to be 0.642 ± 0.013 pj. 

In HQET, the heavy flavor and spin symmetries imply relations among the form factors 
and reduce their number to one when the decay involves only heavy quarks. For heavy-to- 
light transitions, two form factors are needed to describe the hadronic current. In this Letter, 
we follow Ref. |4|, in which the c quark is treated as heavy and the s quark as light. Two 
independent form factors fx and /2 are related to the standard form factors in the following 
way F?tf) = -F?{q 2 ) = h{q 2 ) + ^(q 2 ) and F 2 v (q 2 ) = -F 2 A (q>) = ^/ 2 (<? 2 ). In 
general, / 2 is expected to be negative and smaller in magnitude than f\. If the s quark is 
treated as heavy, / 2 is zero. 

In order to extract the form factor ratio R = /2//1 from a fit to the decay rate, Ts, 
an assumption must be made about the q 2 dependence of the form factors. The model 
of Korner and Kramer (KK) 4| uses the dipole form f(q 2 ) = ^J^h^s ^ (1 - if^) 2 f° r 

' 9 / pole' pole 

both form factors, where the pole mass is taken from the naive pole dominance model: 
M pole =m D , =2.11 GeV/c 2 . 

In this Letter, we perform, for the first time, a simultaneous fit for the form factor ratio 
and pole mass in the decay A+ — > Ae + z/ e , and we make a first search for CP- violation in 
this decay. The data sample used in this study was collected with the CLEO II jfjj and 
upgraded CLEO II. V 0] detectors operating at the Cornell Electron Storage Ring (CESR). 
The integrated luminosity consists of 13.7 fb _1 taken at and just below the T(AS) resonance, 
corresponding to approximately 18 x 10 6 e + e~ — *> cc events. Throughout this paper charge 
conjugate states are implicitly included, unless otherwise indicated, and the symbol e is used 
to denote an electron or positron. 

The analysis is an extension of the technique described in [§, 0|. The decay A+ — ► 
Ae + z/ e is reconstructed by detecting a Ae + pair with invariant mass in the range m^ e + < 
m\ c . The positron is required to come from the region of the event vertex. To reduce the 
background from B decays, we require _R 2 = H2/H0 > 0.2, where Hi are Fox- Wolfram event 
shape variables Positrons are identified using a likelihood function, which incorporates 
information from the calorimeter and dE/dx systems. The minimum allowed momentum 
for positron candidates is 0.7 GeV/c, as the positron fake rates are much higher at lower 
momentum. Positrons must be detected in the region: |cos#| < 0.7, where 6 is the angle 
between the positron momentum and the beam line. The A is reconstructed in the decay 
mode A — > pir~ . The A baryon is long-lived (cr = 7.89 cm), accordingly, the A vertex 
is required to be greater than 5 mm from the primary vertex in the r — <p plane, but the 
A momentum must extrapolate to the primary vertex. The dE/dx measurement of the 
proton is required to be consistent with the expected value. Combinations that satisfy 
interpretation as a Kg are rejected. The magnitude of the A momentum is required to be 
greater than 0.8 GeV/c in order to reduce combinatorial background. These A candidates 
are then combined with Right Sign (RS) tracks consistent with positrons, and the sum of 
the A and e + momenta is required to be greater than 1.4 GeV/c in order to reduce the 
background from B decays. 

The above selection criteria permit the isolation of signal events with low background. 
The number of events passing the selection is 4060, of which 123 ± 12 are consistent with 
fake A background, 338 ± 67 with S c — > Ee + i> feed through and 398 ± 58 with e fake 
background. The sidebands of the pn~ invariant mass distribution are used to estimate the 
fake A background. The background from S c — > Ee + u decays is estimated using the result 
of a previous CLEO analysis pj|. 
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The normalization and momentum spectrum of the e fake background is estimated using 
the Wrong Sign (WS) h + A data sample (no charge conjugation is implied), where h + and A 
satisfy all analysis selection criteria. The h + tracks in this sample are mostly fakes as there 
are few processes contributing e + A pairs after the selection criteria are applied. If no particle 
identification is used for h, the probability to find a (h + A or h~A) WS or (h + A or h~A) RS 
pair is approximately equal because the net charge of the event is zero. The equality is not 
perfect due to (a) baryon conservation: a A is more likely to be produced with an antiproton 
in WS rather than RS combinations, and (b) associated strangeness production: there is a 
higher fraction of kaons in RS rather than WS combinations. When the particle identification 
requirements for h are applied, the importance of these correlations is magnified by the high 
e fake rates of antiprotons and kaons. Therefore, antiprotons are excluded by using only one 
WS charge conjugate state (h + A), and the momentum region where the e fake rate from 
kaons is high is excluded by requiring \p e \ > 0.7 GeV/c. Differences that remain between the 
momentum spectra and particle species of hadronic tracks h + A and h + A are second order 
and are accounted for by a systematic uncertainty. 

Calculating kinematic variables requires knowledge of the A+ momentum, which is un- 
known due to the undetected neutrino. The direction of the A+ is approximated using the 
information provided by the thrust axis of the event and the kinematic constraints of the 
decay. The magnitude of the A+ momentum is obtained as a weighted average of the roots 
of the quadratic equation p A 2 c = (p*A + Pe + Pu) 2 - The weights are assigned based on the 
measured fragmentation function of A+. After the A+ momentum is estimated, the four 
kinematic variables are easily obtained. The kinematic variables t = q 2 /q 2 aax , cos9\, cos 9w 
and x achieve resolutions of 0.2, 0.3, 0.2 and 45°, respectively. 

A four-dimensional maximum likelihood fit in a manner similar to Ref. 0| is performed 
in the space of t, cos^a, cos 8w and \. The technique makes possible a multidimensional 
fit to variables modified by experimental acceptance and resolution taking into account 
correlations among the variables. We have performed two types of fit. The first fit is 
unbinned in all four dimensions. The second fit is unbinned in cos#a, cosQw and x, and 
binned in t. While both fits produce consistent results, the second fit is used for the main 
results. 

The signal probability density function for the likelihood function is estimated at each 
data point using signal Monte Carlo (MC) events, generated according to the HQET consis- 
tent KK model with a GEANT based simulation |l3j], by sampling the MC distribution at 
the reconstructed level in a search volume around the data point. The volume size is chosen 
so that the systematic effect from finite search volumes is small and the required number 
of MC events is not prohibitively high. The background probability density functions are 
modeled similarly using samples of events for each background component. The e fake back- 
ground is modeled using a sample of events collected for that purpose from the data. Feed 
through background from S c — > He + z/ is modeled by the MC sample generated according 
to the HQET consistent KK model. Fake A background is modeled using the data events 
in the sidebands of the pir~ invariant mass distribution. For the binned part of the fit, the 
above distributions are projected onto t and binned. The background normalizations are 
fixed in the fits to the measured values. 

Using the above method, a simultaneous fit for the form factor ratio and the pole mass 
is made. We find R = -0.31 ± 0.05(stat) and M po i e = (2.21 ± 0.08(stat)) GeV/c 2 . This is 
the main result of the analysis. Figures [T] and El show the t, cos9\, cos 6w and x projections 
for the data and fit. 
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FIG. 1: Projections of the data (points with statistical error bars) and the fit (solid histogram) 
onto t, cos 6 a, cosOw and x- The dashed lines show the sum of the background distributions. 




FIG. 2: Projections of the data (points with statistical error bars) and the fit (solid histogram) 
onto cos#a, cosOw and x f° r two t regions. The plots labeled (a), (b) and (c) are for t < 0.5; (d), 
(e) and (f) are for t > 0.5. The dashed lines show the sum of the background distributions. 

We have considered the following sources of systematic uncertainty and give our estimate 
of their magnitude in parentheses for R and M po i e , respectively The uncertainty associated 
with the size of the search volume is measured from a statistical experiment in which a set of 
mock data samples, including signal and all background components, was fit in the same way 
as the data (0.006, 0.048). The uncertainty due to the limited size of the signal MC sample 
is estimated by dividing the sample into four independent equal subsamples and repeating 
the fit (0.007, 0.012). The uncertainty due to background normalizations is determined by 
varying the estimated number of background events by one standard deviation separately 
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for each type of background (0.023, 0.024). The uncertainty associated with the modeling of 
the background shapes, including uncertainties originating from the modeling of the e fake 
background, and the unknown form factor ratio and M po i e for the decays S c — > Se + z/, is 
estimated by varying these shapes or by using alternative background samples (0.024, 0.049). 
The uncertainty due to the small background contribution of random Ae + pairs from the 
continuum (e + e~ — > qq) and T(45) — > BB events, which are not modeled in the fit, is 
obtained from a generic MC sample and is estimated by repeating the fit with and without 
this background (0.013, 0.038). The modes — > AXe + z/, where X represents additional 

decay products, have never been observed. The current upper limit is ^§f^t~^^ < 0.15, 



where X ^ 0, at 90% confidence level |9j. The uncertainty due to the presence of these 
modes is estimated from a series of fits, each having an additional background component 
with floating normalization to represent a A+ — > AXe + v mode. The uncertainty due to the 
possible presence of A+ — > AXe + v is assigned as the largest deviation from the main result 
found in these fits (0.022, 0.091). The uncertainty associated with the A+ fragmentation 
function is estimated by varying this function (0.003, 0.002). The uncertainty related to MC 
modeling of the reconstruction efficiency of slow pions produced in A decays is obtained by 
varying this efficiency according to our understanding of the CLEO detector (0.004, 0.003). 
For approximately 10% of the data, all from runs near the end of the CLEO II. V data 
taking period, a discrepancy between data and MC was found for t > 0.8. The effect of the 
discrepancy was determined by binning the data in t and performing a series of fits with a 
variable range of t excluded. The size of the systematic uncertainty is conservatively taken to 
be the largest difference between the results of these fits and the main result (0.004, 0.072). 
The simulation used to obtain the main result does not include final state radiation. The 
systematic uncertainty due to this is determined as the difference between the result of a fit 
where electroweak radiative corrections have been modeled and the main result (0.005, 
0.009). 

Adding all sources of systematic uncertainty in quadrature, the final result is R = — 0.31± 
0.05(stat)±0.04(syst) and M pole = (2.21±0.08(stat)±0.14(syst)) GeV/c 2 , the latter value be- 
ing consistent with vector dominance. We also find R = — 0.35±0.04(stat) ±0.04(syst) from 
a fit with Mpoio = m D * . Using the values of R and M pole obtained in the simultaneous fit and 
the KK model, the mean value of the decay asymmetry parameter of A+ — ► Ae + u e [3] aver- 
aged over charge conjugate states is calculated to be oa c = —0.86 ± 0.03(stat) ± 0.02(syst), 
for (q 2 ) = 0.67 (GeV/c 2 ) 2 . 

In the Standard Model, CP-violation is expected to be small in semileptonic decays and 
absent in the decay A+ — > Ae + v e . If CP is conserved, the following relation is satisfied 

c/r(A+ -> Ae+z/ e ) dT(A~ -> Ae~V e ) 



dq 2 d cos 9\d cos 9\ydx dq 2 d cos 9\d cos 0wd(— x) 

Following 0| and by extension, a CP- violating asymmetry of the A+ is defined as Aa c = 
("Ac+ a A c ) ^ j n ^ e model, the angular distributions for Aj~ and A are governed by R 

and M pole . Using the values of R and M pole obtained in a simultaneous fit to each charge 
conjugate state separately and the KK model, we calculate «a c «a = —0.561 ± 0.026(stat) 
and cn^ ot^ = —0.544 ± 0.024(stat). From the following expression for the CP asymmetry: 

Ac A . Ac A = .4a + Aa, which is valid to first order in A\ and Aa , we obtain ^4a = 

L A c A 

0.00 ± 0.03(stat) ± 0.01 (syst) ± 0.02, where in the systematic uncertainty the correlations 
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among the systematic uncertainties for the charge conjugate states are taken into account 
and the third error is from the uncertainty in Aa- 

In conclusion, using a four- dimensional maximum likelihood fit, the angular distributions 
of A+ — > Ae + u e have been studied. The form factor ratio R = /2//1 and M po i e are found to 
be -0.31±0.05(stat)±0.04(syst) and (2.21 ± 0.08(stat) ± 0.14(syst)) GeV/c 2 , respectively. 
This is the most precise measurement of R, and it demonstrates that f2 is non zero with a 
combined statistical and systematic significance exceeding 4a. This is also the first measure- 
ment of Mp i e in a charm baryon semileptonic decay. Our measurement is consistent with 
vector dominance. These results correspond to a^ c = —0.86 ± 0.03(stat) ± 0.02(syst), for 
(q 2 ) = 0.67 (GeV/c 2 ) 2 . Comparing the angular distributions for A+ and A c , no evidence for 
CP-violation is found: A Ac = 0.00 ± 0.03(stat) ± O.Ol(syst) ± 0.02, where the third error is 
from the uncertainty in the world average of the CP- violating parameter, Aa, for A — > pn~ . 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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